Mon. Not. R. Astron. Soc. 000, 000-000 (2007) Printed 5 February 2008 



(MN WF$t style file v2.2) 



Imaging and Spectroscopy of Ultra Steep Spectrum Radio 
Sources *. 

Carlos G. Bornancini 1 ' 2 , Carlos De Breuck 3 , Wim de Vries 4 ' 5 , Steve Croft 4 ' 6 , Wil 
van Breugel 4 ' 6 , Huub Rottgering 7 and Dante Minniti 8 

1 Grupo de Investigaciones en Astronomia Teorica y Experimental, IATE, Observatorio Astronomico, Universidad National de Cordoba, 
Laprida 854, X5000BGR, Cordoba, Argentina. 

2 Secretaria de Ciencia y Tecnica de la Universidad National de Cordoba. 

3 European Southern Observatory, Karl Schwarzschild Strafie 2, D-85748 Garching, Germany. 

4 Institute of Geophysics and Planetary Physics, Lawrence Livermore National Laboratory L-413, 7000 East Avenue, Livermore, CA 94550, USA 

5 University of California, Davis, 1 Shields Avenue, Davis, CA 95616, USA 

6 University of California, Merced, P. O. Box 2039, Merced, CA 95344, USA 

7 Leiden Observatory, PO Box 9513, 2300 RA, Leiden, The Netherlands 

8 Pontificia Universidad Catolica de Chile, Departamento de Astronomia y Astrofisica, Casilla 306, Santiago 22, Chile 
5 February 2008 

ABSTRACT 

We present a sample of 40 Ultra Steep Spectrum (USS, a < —1.3, S v oc v a ) radio 
sources selected from the Wcstcrbork in the Southern Hemisphere (WISH) catalog. 
The USS sources have been imaged in i^-band at the Cerro Tololo Inter-American 
Observatory (CTIO) and with the Very Large Telescope at Cerro Paranal. We also 
present VLT, Keck and Willian Herschel Telescope(WHT) optical spectroscopy of 14 
targets selection from 4 different USS samples. For 12 sources, we have been able to 
determine the redshifts, including 4 new radio galaxies at z > 3. We find that most 
of our USS sources have predominantly small (<6") radio sizes and faint magnitudes 
(K>18). The mean if— band counterpart magnitude is i\~=18.6. The expected redshift 
distribution estimated using the Hubble K — z diagram has a mean of z exp ^2.l3, which 
is higher than the predicted redshift obtained for the SUMSS-NVSS sample and the 
expected redshift obtained in the 6C** survey. The compact USS sample analyzed here 
may contain a higher fraction of galaxies which are high redshift and/or are heavily 
obscured by dust. Using the 74, 352 and 1400 MHz flux densities of a sub-sample, we 
construct a radio colour-colour diagram. We find that all but one of our USS sources 
have a strong tendency to flatten below 352 MHz. We also find that the highest redshift 
source from this paper (at z=3.84) does not show evidence for spectral flattening down 
to 151 MHz. This suggests that very low frequency selected USS samples will likely be 
more efficient to find high redshift galaxies. 

Key words: surveys - radio continuum: general - radio continuum: galaxies - galax- 
ies: high-redshift 



1 INTRODUCTION 



* Based on observations obtained at Cerro Tololo Inter-American 
Observatory, a division of the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. Based on observations with the 
Very Large Telescope, obtained at the European Southern Obser- 
vatory in Chile under Proposals 65.O-0125(A), 66.A.0006(B) and 
69.A-0337(A). 



Distant powerful radio sources represent excellent tar- 
gets to study evolutionary processes related to mas- 
sive galaxies and their associated surrounding struc- 
high 



tures at high r edshifts (e.g. |Pe Breuck et afT l2002bl; 
Miley et al.1 12004 lOverzier et al.1 [20061 : IVillar-Martfn et al.l 
20031 : iNesvadba et alj|2006r i. The most effective method for 
finding high-redshift radio galaxies is by selecting for Ultra 
Steep-Spectrum radio sources (USS, a < —1, S v oc v a ). 
A possible explanation of this empirical method is a K— 
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correction induced by a curvature of the radio spectra. 
High redshift sources with concave radio spectra will have 
steeper spectral indice s compared to those at low redshift 
llKrorik fc Chenl Il99ll ; ICarilli et al] Il999l ; iGopal-Krishnal 
1988). Comparing the extremely steep spectral index sources 
assoc iated with galaxies r esiding closest to the cluster cen- 
tres, iKlamer et al.l (|2006l ) found an explanation that sug- 
gests that the steeper spectra can also be explained by 
pressure-confined radio lobes which have slow adiabatic ex- 
pansion losses in high density environments. High-redshift 
radio galaxies are extremely luminous, and spatially ex- 
tended, compared to normal galaxies at similar redshifts. 
The tight K — z relation of radio galaxies in the Hubble K 
diagram suggest that these object s are amongst the most 
massive systems a t each redshift toe Breuck et aDl2 002b; 
Jarvis et al. l200ll; IWillott et all 120031 ; lEales fc Rawlingsl 
19931 ; Lacy et alTl200d ). Using the K — z diagram and 



galaxy evolution models, iRocca-Volmerange et alJ (|2004l ) 
found that the typical hosts of radio galaxies correspond 
to the most massive elliptical galaxies with baryon masses 
M ~ 1O 12 M . 

There is increasing new evidence that d istant radio 
;alaxi es represent massive forming systems. iMilev et al.l 
2006) found more than 10 individual clumpy features, pos- 
sible satellite galaxies in the process of merging, in a re- 
gion of 50x40 kpc around the radio galaxy MRC 1138-262 
at z = 2.2. From a large VLT program, searching Lya 
emitters around a sam ple of 2 < 2 < 5.2 radio galaxies, 
IVenemans et al.1 l|2006l ) found that at least six of eight fields 
studied are overdense in Lya emitters by a factor of 3-5 
compared to the field. Some radio-galaxies, selected with 
the USS criterion, are surrounded by giant Lya haloes. The 
size of these gas structures (100-200 kpc) are similar to the 
observed size of cD haloes in galaxies in ne arby clusters 
l|Reuland et al.ll2003al ; Ivan Breugel et aLlbOOrj ). From a spa- 
tial cross-correlation analysis betw een USS at 0.5 < z < 1.5 
and surrounding K < 20 galaxies, iBornancini et all (120061 ) 
found a comoving correlation length comparable to those ob- 
tained for clusters of galaxies with masses in the range M ~ 
W 14 M Q in a cosmological N-body simulation. In this paper 
we present if-band observations, high-resolution radio maps 
and optical spectroscopy for a subsample of southern Ultra 
Steep Spectrum Radio Sources selected f rom the Wester- 
fa ork in the Southern Hemisphere (WISH, |Pe Breuck et al] 
(2002a)) 352 MHz catalog. We als o present optical spec - 
troscopy of sources selected from the lDe Breuck et al.l (|2000h 
and WISH USS samples, yielding 12 new redshifts, includ- 
ing four new radio galaxies at z>3 These provide excellent 
targets for follow-up observations in this interesting redshift 
range. 

The structure of this paper is organized as follows: Sec- 
tion 2 describes the sample analyzed. We describe the near- 
IR observations and data reductions in Section 3. Section 4 
presents the near-IR identification and source extraction. We 
describe optical spectroscopy observations in Section 5 and 
discuss our results in Section 6. Finally, Section 7 presents 
the main conclusions. 

In this work we assume a standard ACDM model Uni- 
verse with cosmological parameters, S7m=0.3, S7a=0.7 and 
a Hubble constant of H =100 Km s -1 Mpc"\ 



2 SAMPLE DEFINITION 

The USS sample was selected from the 352 MHz Westerbork 
In the Southern Hemisphere catalogue (|De Breuck et al.l 
2002a). We used a selection criterion based on the ra- 
dio properties of the sources, steep radio spectral index 
( a 352° < —1-3) and small angular size (< 30", as measured 
on the VLA maps). All of the first and most of the second 
priority targets have been observed. A detailed description 
of the radio observ a tions a nd data reduction can be found 
in |Pe Breuck et all (|2002al ). Most of the targets selected for 
optical spectroscopy ar e selected from the USS samples of 
iDe Breuck et all (|200Cj ). 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Radio Imaging 

3.1.1 VLA 

In order to obtain high resolution radio maps and accurate 
positions, we used the Very Large Array (VLA) of the Na- 
tional Radio Astronomy Observatory in the hybrid BnA- 
configuration at 1.4 GHz. The images were obtained on 1999 
October 16 and 20 and consisted of L-band short subscans of 
3 minutes each. We used the standard data reduction recipes 
in the Astronomical Image Processing Software (jAX'PS ) 
package, including self-calibration for phase corrections. 



3.2 K-band Imaging 

3.2.1 CTIO 

The USS fields were observed during two runs in 2000 
March and 2001 January using the Ohio State InfraRed Im- 
ager/Spectrometer (OSIRIS) imager on the 4— meter V.M. 
Blanco telescope at the Cerro Tololo Inter- American Ob- 
servatory (CTIO). OSIRIS is a 0.9-2.4 ^m camera with a 
1024x1024 HAWAII HgCdTe CCD array. We used a focal ra- 
tio of f/7 which resulted in a respective pixel size of CK' 161. 
Individual frames were obtained as a co-addition of a num- 
ber of single exposures with different integration times, in 
this way one effectively removes row/column related defects. 
For example, for sources with total exposures times of 1920s, 
we use 12 single exposures of 10s observed with a a 16— point 
dithering pattern. This results in 16 x 10 x 12 = 1920s = 32 
minutes on source. We used similar procedure for the other 
images, doing small changes in the number of points of the 
dithering pattern and the exposure times for single obser- 
vations During the observations, tip-tilt mirror corrections 
were made. This resulted in typical FWHMs on the individ- 
ual pointing frames (12 x 10s co-added) of 0'.'5 to 0'.'7. The 
mean was around tf!6 for the nights 2000 March 20 to 22, 
and around 0'.'7 on the 2001 January run. 

A subsample of USS sources was observed in 1999 July 
with the Cerro Tololo IR imager (CIRIM) camera at the 
CTIO 4-meter V.M. Blanco telescope. The CCD detector 
is a 256x256 HgCdTe NICMOS 3 array with a pixel scale 
of 0'.'414 using a focal ratio of f/7. 5. For these observations 
we obtained a mean FWHM of l'.'l. 
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3.2.2 VLT 

Seven sources undetected in the CTIO images were observed 
with the Infrared Spec trometer And Array Camera (ISAAC; 
iMoorwood et al.lll99sh on the Very Large Telescope (VLT) 
at Paranal, (Chile) between 2002 April and September. A 
K a filter (2.0-2.3 /jm) was used. The pixel scale was 0'.'I48 
per pixel . All the images were obtained under optimum see- 
ing conditions with FWHM in the range Of! 7 and l'.'O, and 
a number of standard stars were observed during the same 
nights. In Table 1 we present the log of the near-IR observa- 
tions, with the observation dates, telescope/instrument used 
and the total exposure times for all the USS sample. 

3.2.3 Data Reduction 

Data were reduced within irafQ, using the standard DIMSUM 
l|Stanford et al.lll995l ) near-IR reduction package, including 
dark subtraction, flat fielding, sky-substraction, bad pixel 
masking, bright object masking, registration and summing. 

VLT images were reduced with ISAAC pipeline recipes 
5.4.23 using the command line utility EsoRex, which includes 
standard routines, including dark correction, flat field cali- 
bration, bad pixel cleaning, image correlation and recombi- 
nation. 

3.2.4 Astrometry 

The astrometric calibration of the images was performed 
using the WCStools packagfjl (|Minkll20od ). using the Two 
Micron All Sky Surv ey (2MASS) Cata logue of Point Sources 
a reference catalog JCutri et alJl20o3 ) and the USNO-A2.0 
Catalogue l|Monet et alj 19981 ) for fields with a few number 
of 2MASS sources. We estimate the uncertainty in the rela- 
tive astrometry to be ~0'.'4-0'.'6 . 

3.2.5 Photometry 

The photometry for the CTIO observations was cali- 
brated using 2MASS point sources detected in these fields. 
For VLT runs, magnitudes we re calibrated using NIC- 
MOS near- infrared standards l|Persson et all [l998). All 
standard star mag nitudes were obtained using SExtractor 
l|Bertin fc Arnoutall996h with an aperture diameter of 6", 
which proved to be adequate by monitoring the growth curve 
of all the measured stars. For the observations made with 
OSIRIS and CTRIM at CTIO the A"-band zero point was 
22.64±0.06 and 22.01±0.02 (for 1 count/second, integrated 
over the source), respectively. For the ISAAC/ VLT obser- 
vations we derive the zero-point based on observations of 
the near-I R photometric s tanda rd star S301-D taken from 
the list of iPersson et al. I (| 19981 ). The zero-point obtained 
was 24.18±0.01. The zero-points do not include the airmass 
term, because the airmass dependence in 7^-band is small 
and we have observed all our objects with airmasses < 1.3, 
with a mean of 1.06. 

1 Image Reduction and Analysis Facility (IRAF), a software sys- 
tem distributed by the National Optical Astronomy Observatories 
(NOAO) 

2 Available at ftp://ftp.eso.org/pub/dfs/pipclincs/isaac/ 

3 Available at ftp://cfa-ftp.harvard.edu/pub/gsc/WCSTools/ 



Because of the large uncertainties in the photometry, 
we did not correct the A'-band magnitudes for Galactic 
extinction. We have used the NED Databas^B, which are 
based on the E(B — V) values from the extinction maps of 
ISchlegel et alj (|l99cf ). We found that most of the sources 
have values from 0.01 to 0.1, which are negligible compared 
to magnitude errors. 



4 SOURCE DETECTION AND NEAR-IR 
IDENTIFICATION 

We have identified the near-IR counterparts to the radio 
sources by overlaying the radio contour maps on to the K— 
band images using the AXPS task KNTR. 

We measured the magnitudes of these identifications 
using SExtractor. The parameters were set such that, to be 
detected, an object must have a flux in excess of 1-1.5 times 
the local background noise level over at least N connected 
pixels, according to the seeing conditions and image quali- 
ties. In order to improve the detections of faint sources, we 
smoothed the Tf-band images using a circular Gaussian of 
FWHM=l-2 pixels. In some cases, we find that it is possi- 
ble that the identification does not fall at the midpoint of 
the radio lobes, for radio sources with double or multiple 
components. 

We measured magnitudes for all radio source identi- 
fications using SExtractor and circular apertures with di- 
ameters of 2'.'0, 4'.'0 and 8'.'0. We also obtained MAG_BEST 
magni tudes in order to facilitate co mparison with the liter- 
ature (|Bornancini et all |2004| . I2OO6I ). This is usually equal 
to MAG_AUT0, but if the contribution of other source exceeds 
10%, it is MAG_IS0C0R. The uncertainties quoted represent 
our best estimates, which includes the zero-point uncer- 
tainty (See Table 2). 

Table 2 lists our USS sample, the columns are: IAU 
J2000 designation, spectral index obtained between 352 
MHz and 1.4 GHz, largest angular size, as measured from ra- 
dio maps taken with the VLA, aperture magnitudes in 2'.'0, 
4'.'0 and 8'.'0, SExtractor MAG_BEST magnitudes, J2000 posi- 
tion coordinates of the radio source and the K-banA iden- 
tification counterpart, the half-light radius and CLASS_STAR 
values of the A"-band identifications. SExtractor CLASS_STAR 
parameter provides an indication of the likelihood of an ob- 
ject being a galaxy or a star. In the ideal cases a galaxy has 
CLASS_STAR= 0.0 and a point source has CLASS_STAR= 1.0. 
We find that all the USS sample analyzed in this work 
are consistent with extended sources (galaxies) rather than 
point sources (stars or quasars). 

We also considered studying the environments of these 
USS radio galaxy candidates. The environments of twenty 
USS sources observed w ith CTIO/OSIRIS were analyzed in 
iBornancini et all l|2004l ). In those images, only 7 radio galax- 
ies were iden tified, finding a statistica lly significant signal of 
clustering. In lBornancini et al.l ()2006h . a similar analysis was 
performed on a wider field and with deeper data, finding a 
stronger result. We therefore do not analyze the environ- 
ments of the USS sources in the relatively shallow K— band 
images from this paper. 

4 http://nedwww.ipac.caltech.edu/ - the NASA-IPAC Extra- 
galactic Database 
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5 OPTICAL SPECTROSCOPY 

During the past 6 years, we have obtained optical spec- 
troscopy of 14 USS sources. The targets were mostly se- 
lected from 4 different USS samples: the 3 sub-samples of De 
Breuck et al. (2000) and the WISH-NVSS sample described 
in this paper. From these samples, we mostly selected the 
final targets based on the RA range when telescope time was 
available. 

5.1 VLT 

Optical spectroscopy of a subsample of USS sources were 
carried with FORS2 installed on the VLT Unit Telescope 

4 Yepun between 2000 October a nd 2001 May. We used a 
similar data reduction strategy to |Pe Breuck et all |2006). 
We removed cosmic rays using the IRAF task szap, af- 
ter bias and fiatfield correction. We used the IRAF task 
background in order to remove sky variations and we ex- 
tracted the spectra using a width appropriate to contain all 
flux in the extended emission lines. Finally we calibrated the 
one dimensional spectra in wavelength and flux. We used 
th e procedures des c ribed in iRottgering et alj (|l997t ) and 
in |Pe Breuck et al] l|200rj ) to determine the central wave- 
length, total line flux, deconvolved widths, and rest-frame 
equivalent widths (with their related errors) of all the emis- 
sion and absorption lines detected in our spectra. For each 
object, we quote a single redshift, which is an average of 
all the lines, weighted by the uncertainties due to the line 
fitting and wavelength calibration. In Table 4 we show the 
measured parameters of the emission/absorption lines in the 
spectra. 

5.2 Keck 

Longslit observations were made on UT 2004 January 19 - 
20, using the Lo w Resolution Imaging Spectrometer (LRIS; 
lOke et all 1 19951 ) on Keck I, with a slit width of l'.'5, un- 
der photometric conditions and / .'8 seeing. The spectro- 
scopic setup is described in Table 3. Data reduction was 
performed in IRAF; data were bias-subtracted, response- 
corrected, cosmic-ray-zapped using Pieter van Dokkum's 
LACosmic routine 0, sky-subtracted, and spectra extracted 
in a l'.'5 aperture. The resultant 1-D spectra were wave- 
length and flux calibrated using appropriate arcs and stan- 
dard stars. 

5.3 WHT 

On UT 2006 December 14 and 15, we obta ined 3 spec- 
tra u sing the ISIS double beam spectrograph |Carter et all 
1 1994 ) at the 4.2 m William Herschel telescope (WHT) at 
Observatorio Roque de los Muchachos on the Canary is- 
land of La Palma. Conditions were photometric with 0'.'5 
and 1'.'5 seeing during the first and second night, respec- 
tively. We used a dichroic splitting the light at 5300 A and 
a l'/5 wide slit. The grisms used were the R300B in the 
blue arm and the R316R in the red arm. None of the 3 ob- 
jects were detected in the blue arm. In the red arm, only 

5 http:/ /www. astro. yale.edu/dokkum/lacosmic/ 



WN J0912— 1655 was clearly detected in 5400 s of observing 
time. WN J0604-2015 and WN J0610+6611 were not de- 
tected in 5400 s and 7200 s, respectively. We extracted the 
red spectrum of WN J0912— 1655 with an aperture width of 
2'.'4. 

5.4 Results 

We have determined redshifts of 12 out of 14 targets ob- 
served. The two targets without redshifts were undetected 
in 1.5 to 2 hours on a 4 m telescope, and need long integra- 
tion on 8-10 m class telescopes. Apart from TN J2009-3040, 
which is a 2=3.158 quasar, all 11 targets are radio galaxies 
with redshifts from z=0.550 to z=3.837. Five of the radio 
galaxies are at z>2, including three at z>3. Table 4 lists 
the emission line parameters of all targets. The 10 new red- 
shifts from the De Breuck et al. (2000) USS sample raises 
the number of spectroscopic observations from that sample 
from 46 to 56, and the number of known redshifts from 34 
to 44. 



6 DISCUSSION 

6.1 Correlations 

In this work we have excluded large radio sources because 
those are most likely foreground objects and not at high 
redshift. In Figure 1 we present the correlation between 
largest angular size (LAS) and the K— magnitudes obtained 
within an aperture of 4". As can be seen most of the sources 
have small radio sizes with LAS < 8". The vertical dashed 
line represent the expected K— band magnitude of a z ~ 3 
source, indicating that high redshift sources have compact 
radio morphologies. Several authors have applied the small 
angular- size criterion in order to find high redshift radio 
galax ies (|Blundell et al.lll998l ; rJarvis et al.ll2004l ; ICruz et all 
[2006|). There are a few cases of radio galaxies at high red- 
shifts with large ang ular siz e s (4C 23.56 at z=2.483 has 
LAS=53"). However, IPedanil (|2003h found that the intro- 
duction of a strong angular size bias such as LAS< 15", in- 
creases the efficiency in selecting high redshift radio galaxies 
by a factor of two. 

Some radio galaxies have very faint or no A-band 
counterpart presenting large deviati on in the Hubble K — 
z diagram l|De Breuck et all |2006h . and 20% of high- 
redshift radio galaxies fail to show optical emission lines 
in spectroscopy observations using the Ke ck Telescope 
jDe Breuck et all l200ll ; iReuland et all l2003bh . This could 
be an indication that some h igh redshift radio ga laxies con- 
tain great amounts of dust. IWillott et all l|2002l ) found an 
anti-correlation between dust emission (mea sured from 850 
umflux) and the projected linear size, though Fkeuland et all 
|2004h did not confirm this using a larger sample of radio 
galaxies. It is likely that the compact USS sample analyzed 
here may also contain a mix of galaxies either at very high 
redshift, or heavily obscured by dust. 

6.2 Expected redshift distribution 

We have estimated the redshift of the USS sample us- 
ing linear regression and the Hubble K — z diagram of 
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Figure 1. Radio largest angular size vs K— magnitudes (4" aper- 
ture). Dashed line represent the expected K— band magnitude of 
a z ~ 3 source. Most of the sources have small radio sizes as 
indicated by the horizontal dotted line. 




2 3 4 

Expected redshift 



Figure 2. Distribution of expected redshift of USS sources, es- 
timated from the Hubble K — z diagram using 64 kpc metric 
apertures. 



iDe Breuck et al.l l|2002bJ ). Using aperture p hotom etry in 
a 64 kpc metric aperture, IDe Breuck et all l|2002al ) found 
K =4.633 xlogio (z) — 17.266. As we expect most radio galax- 
ies to have 1 < z < 4, we have assumed K^ikpc = K (8"), as 
8" corresponds to ~64kpc at 2=1, which only a small depen- 
dence on redshift z > 1 in our adopted Cosmology. Figure 2 
shows the expected redshift distribution for our USS sample. 
The median predicted redshift is z exv =2.1 3, which is higher 
than the predicted redshift obtained by IDe Breuck et al] 



20041 ) for the SUMSS-NVSS sample (z sumss = 1.75) 
and the expected redshift obtained in the 6C* (z ~ 1.9; 

1.7; ICruz et all 



pe 

Jarvis et al.ll200l] ) and 6C** surveys (z 
20071 ). 



In Figure 4 we show the K — band magnitude distribu- 
tion of USS sources, measured in a 8"diameter aperture. 
We co mpare our results with t hose ob tained in the SUM SS- 
NVSS (|De Breuck et al.ll2004l) . 6C** dCruz et al.ll2007l) cat- 
alogue s and data selected at 74 MHz taken from lJarvis et alj 
l|2004h . The mean K— band magnitudes in an 8" diameter 
aperture of the USS is A=18.5 8, which is one magnit ude 
fainter than those obtained by |Pe Breuck et all (|2004|) m 



the S UMSS-NVSS survey (#=17.57) and by ICruz et all 
|2006l ) in the the 6C** survey (A=17.59). We find a sim- 
ilar result using BEST magnitudes. Our USS sample has 
Abest = 18.74, and from a similar study iBornancini et all 

(2006) found A B est=17.97. The fainter A-band counterparts 
of the WISH-NVSS USS sources could be an indication that 
the host galaxies are located either at higher redshifts (i.e. a 
more efficient selection), or in a very dense, dusty medium. 

Alternatively, the WISH-NVSS USS radio galaxies are 
intrinsically less massive and hence less luminous than the 
ones in previous samples. This could be due to the fainter 
flux limits probed in this sample. However, using Spitzer 
rest-frame near-IR phometry, Seymour et al. 2007 do not 
find a strong dependence of host galaxy mass on 3 GHz radio 
luminisity, so we regard this explanation as unlikely given 
the small difference in radio luminosity between the USS 
samples. 

Alternatively, the WISH-NVSS USS sources may be 
less massive and hence fainted in the observed A— band. 
This may be because WISH-NVSS sources have fainter radio 
fluxes than most previous USS samples. Such a radio power 
dependence i n the K — z diagram has been suggested by var- 
ious authors ([Eales et al.lll997|:lBest et alJll99ct|Jarvis et alj 
l200ll ; |Pe Breuck et al.ll2002bl ; Iwillott et al.ll2002l ). However, 
using rest-frame near-IR imaging of a sample of 69 radio 
galaxies at z > 1 observed with Spitzer, Seymour et al. 

(2007) do not find a strong dependence of host galaxy mass 
on 3 GHz radio luminisity, so we regard this explanation as 
unlikely given the small difference in radio luminosity be- 
tween the USS samples. 

6.3 Radio two-colour diagram and spectral 
curvature 

We use th e CATS database of the Special Astronomy Ob- 
servatory (jVerkhodanov et al.lll997T ) to search for radio mea- 
surements in order to study the radio spectral energy dis- 
tribution at low frequencies. Figure 3 shows the radio two- 
colour diagram, which compares the spectral indices at 74- 
352 MHz and 352-1400 MHz for 12 USS sources of the total 
sample (Filled circles) . Fluxes at 74 MH z were obtained from 
the VLA Low- frequency Sky Survey ijCohen et all I2006T ). 
The line indicates the relation for USS sources whose spectra 
follow a single power law from 74 to 1400 MHz. Open circles 
represent measurements obtained using th e WENSS-NVSS 
sampl e at 325 MHz and 1.4 GHz taken from lPe Breuck et al.l 
(2000). We find a general agreement with those values. Only 
a single source in our sample is consistent with a straight 
spectrum, and it is clear that most of the USS sources have 
substantially flatter spectral indices between 74 and 352 
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Figure 3. Radio two-colour diagram for USS sources in our sam- 
ple (filled circles). Open circles represent measurements obtained 
usin g the WENSS-NVSS sam ple at 325 MHz and 1.4 GHz taken 
from [be Breuc k et al. The line indicates the relation for 

USS sources whose spectra follow a single power law from 74 to 
1400 MHz. 
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Figure 4. K— band magnitude distribution of USS sources, mea- 
sured in a 8"diameter aperture (shaded histogram). Solid and 
dashed lines represent the distribution of K— band magnitude, 
measured in a 8 "diameter aperture for sources dete cted in the 
SUM SS-NVSS jDe Breuck et al.l |2004|) and 6C** llCruz et al.l 
l2007h catalogues. Dotted histogram represent USS sources se- 
lected at 74 MHz taken from ljarvis et alj j2004h . 



MHz than between 352 and 1400 MHz. The median spec- 
tral index obtained at 74-352 MHz is S74 2 = -1.1, while 
the median at high frequencies is S352 = —1.5. The spec- 
tral flattening is in the range 0<A(a)<0.8. Unfortunately, 
we do not have redshifts for any of the 12 sources with 
74 MHz data. However, WN J2007-1316 at z=3.837 is de - 
tected in the Mauritius Southern Sky Survey (|Pandevll2006l ) 
with 5'i5iMHz=3.58±0.57 Jy. This implies a?|?=-1.60±0.19, 
compared to aj^ — — 1-52±0.04. Although the uncertainty 
in the spectral index is larger than for the 74-352 MHz data, 
this result suggests that the highest redshift radio galaxy 
confirmed in our sample does retain a straight radio spec- 
trum through the lo west frequencies observed. As noted by 
iKlamer et all l|2006l ) ~90% of USS sources selected from the 
SUMSS and NVSS catalogues show a straight radio spec- 
tra characterised by a single power law beyond ~1 GHz. 
At frequencies below 100 MHz spectral curvature is much 
more common than at higher frequencies due principally to 
synchrotron self-absorption. This suggests that very low fre- 
quency selected USS samples will be likely to be more effi- 
cient to find high redshift galaxies. 



7 CONCLUSIONS 

We present Jf-band observations obtained with CTIO 
at Cerro Tololo and with VLT at Cerro Paranal, high- 
resolution VLA radio maps for a sample of 28 Ultra Steep 
Spectrum Radio Sources selected from the WISH catalog. 
We also present optical spectroscopy from VLT/FORS2, 
Ke ck/LRIS and WHT/ISIS for 12 USS sources, mainly from 
the lDe Breuck et al.l i|2000T ) samples. 

We find that most of our USS sources have small ra- 
dio sizes and fainter magnitudes. The mean if —band coun- 
terpart magnitude measured in 8" diameter aperture is 
7f=18.58, which is one magnitude fainter than those ob- 
tained on the literature. The expected redshift distribution 
estimated using the Hubble K — z diagram has a mean of 
2cip=2.13, which is higher than the predicted redshift ob- 
tained for the SUMSS-NVSS sample and the expected red- 
shift obtained in the 6C* and 6C** surveys. 

From our spectroscopy, we identify one 2=3.158 quasar, 
and 11 new radio galaxies, 3 of which are at z>3. 

From a radio colour-colour diagram of a subsample of 
12 sources, we find that all but one radio spectra flatten 
significantly below 352 MHz. We also find that the highest 
redshift source from this paper (at 2=3.84) does not show 
evidence for spectral flattening down to 151 MHz. Because 
most low redshift radio galaxies have a flatter spectra below 
352 MHz, selecting USS sources at low frequencies is likely 
to be more efficient selection technique to find high redshift 
galaxies. The forthcoming low frequency observatories such 
as LOFAR, LWA or SKA will be enable to f ind large samples 
of low - frequency selected USS sources fe.g. lRottgering et al.l 
2006) and will therefore open new perspectives on the nature 
and evolution of high redshift radio galaxies. 
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Figure 5. Part of the i_"-band images of the WISH sample from observation of CTIO/OSIRIS, CTIO/CIRIM or VLT/ISAAC, with 
radio contours overlaid. The contours represent the VLA radio emission at 1.4 GHz. The contour scheme is a geometric progression in 
y/2, which implies a factor two change in surface brightness every 2 contours. The first contour level is at 3<7 rm s, where <7 rm s is the rms 
noise measured around the sources, which indicated above each plot. The half-power beamwidth is indicated in the lower left corner of 
the plots. The open cross indicates the K— band identification, as listed in Table 2. See |http://www. eso.org/ ~cbreuck/papers.html| for 
a version with all 28 overlays 



which improved the paper. We thank Adam Stanford and 
Cyril Tasse for their help in taking the Keck and WHT spec- 
tra, respectively. This work was partially supported by the 
Consejo Nacional de Investigaciones Cientfficas y Tecnicas 
(CONICET), the Secretaria de Ciencia y Tecnica de la 
Universidad Nacional de Cordoba and the European Union 
Alfa II Programme, through LENAC, the Latin American- 
European Network for Astrophysics and Cosmology. The 



work of WdV, WvB and SDC was performed under the 
auspices of the U.S. Department of Energy, National Nu- 
clear Security Administration by the University of Califor- 
nia, Lawrence Livermore National Laboratory under con- 
tract No. W-7405-Eng-48. Based on observations made with 
the European Southern Observatory telescopes obtained 
from the ESO/ST-ECF Science Archive Facility. The au- 
thors made use of the database CATS (Verkhodanov et 



8 C. Bornancini et al. 



MP J0340-6507 z=2.289 



WN J0528+6549 z=1.215 




6 ; 

1 4 | 

: 



5000 6000 
Wavelength (Angstrom) 

MP J0601-3926 z=1 .633 

cm] cii] 




5000 6000 
Wavelength (Angstrom) 

WNJ0912-1655 z=0.910 

[oil] 



6000 7000 
Wavelength (Angstrom) 

WN J1 053+5424 z=3.083 



4500 5000 
Wavelength (Angstrom) 

WNJ2007-1316 z=3.837 



CIV Hell 



7000 8000 
Wavelength (Angstrom) 

MPJ2226-7654 z=0.550 



7000 8000 
Wavelength (Angstrom) 









7500 8000 
Wavelength (Angstrom) 

WN J0716+5107 z=1.140 




7000 7500 8000 

Wavelength (Angstrom) 

TN J1049-1258 



6 
4 

I 2 



-4 k 
4500 



5500 6000 6500 

Wavelength (Angstrom) 



TNJ1941-1951 z=2.667 



4000 5000 6000 7000 

Wavelength (Angstrom) 



TN J2009-3040 z=3.158 



5000 6000 7000 

Wavelength (Angstrom) 

MPJ2352-6154 z=1.573 




7000 7500 





S 4 
—i 

„> 2 




5000 6000 
Wavelength (Angstrom) 



Figure 6. Optical spectroscopy of USS sources. Prominent features are indicated with vertical dotted lines showing the predicted 
observed wavelength of the lines at the redshift the source. Atmospheric absorption features are indicated with ©. 
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Table 1. Log of the /\-band observations of the USS WISH sample. 
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Table 2. Results of the if-band and radio observations of the WISH-NVSS USS sample 
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Table 3. Log of the Keck I optical spectroscopy. 



Target 


Dichroic 


Blue grism 


Blue pix. scale 


Resolution 


Red grating 


Red pix. scale 


Resolution 


Exp. time 


PA 


WN J0528+6549 


D680 


400/3400 


1.04 A / pix 


9.5A 


400/8500 


1.86A / pix 


8.lA 


1 X 1800 s 


177.00 


WN J0716+5107 


D680 


400/3400 


1.04 A / pix 


9.5A 


400/8500 


1.86A / pix 


8.lA 


2 x 1800 s 


201.00 


WN J1053+5424 


D560 


300/5000 


1.43 A / pix 


12A 


400/8500 


1.86A / pix 


8.lA 


4 x 1800 s 


254.06 



Table 4. Emission line measurements 



Source 


z 


Line 


Aobs 


Flux 


A"FWHM 










A 


10~ 16 crgs _1 cm -2 


km s _1 


A 


TV/TP Ifi'iAn fi^fl7 


9 90Q_i_n nn^ 
z.zoyzcu.uuo 


Lya 


AC\f\C\~\-A 
'IUUUZl 1 


n ziri-t-fi 

U.^tUICU.UO 




,>OG 






U IV Alo4y 


c:i nn4-07 
OlUUztZ l 


u. lOztu.UO 


9Qnri-i- 1 son 


on_i_ 1 n 






He II A1640 


6275±10 


0.26±0.05 


1400±500 


>40 






C III] A1909 


6270±10 


0.25±0.05 


1600±1100 


>130 


vv duozoTUO^y 


1 91 7-i-n nn^ 
i.zi / znu.uuo 


m TTl \Q797 

[\J 11J AO ( Z / 


OZUOZL ± t 


n HQ-l-n 09 






MP J0601-3926 


1.633±0.001 


C IV A1549 


4083±27 


0.30±0.15 


2300±1600 


>4 






C III] A1909 


5020±4 


0.6±0.1 


1200±500 


28±6 






P TTl \9^9£ 
11J AZoZO 


OloUzC 1 




OOUztoDU 




WN J0716+5107 


1.1401±0.0002 


[0 II] A3727 


7976±1 


0.06±0.01 


<500 


8±1 


WN J0912-1655 


0.9102±0.0002 


[O II] A3727 


7119±1 


0.10±0.02 


<100 


11±2 


TN J1049-1258 


3.697±0.004 


Lya 


5712±5 


0.6±0.1 


1600±500 


>63 


W1V T1 fi^-i-^AOA 
VV IN J lUOOTi)^^ 


o.uooznu.uui 


Lya 




u.uuznu.ui 






TN J1941-1951 


2.667±0.001 


Lya 


4460±1 


3.7±0.4 


800±200 


330±70 






C IV A1549 


5680±2 


0.9±0.1 


850±200 


64±20 






He II A1640 


6010±3 


0.4±0.1 


850±300 


23±6 






C III] A1909 


7000±16 


0.2±0.1 


1000±700 


10±4 


WN J2007-1316 


3.837±0.001 


Lya 


5883±3 


1.8±0.2 


2400±300 


77±13 






C IV A1549 


7500±20 


0.5±0.1 


3100±1700 


26±5 






He II A1640 


7933±2 


0.2±0.1 


450±200 


9±2 


TN J2009-3040 


3.158±0.001 


Lya 


5056±1 


15.8±1.6 


2300±500t 


65±7 






C IV A1549 


6444±1 


2.0±2.0 


2300±500t 


10±1 






He II A1640 


6815±5 


0.65±0.1 


1380±500 


4±1 






C III] A1909 


7920±5 


3.0±0.4 


2900±600t 


14±2 


MP J2226-7654 


0.5500±0.0002 


[O II] A3727 


5776±1 


2.8±0.3 


400±180 


25±3 






Call K 


6098±1 


-2.3±0.3 


1200±300 


-12±2 






Call H 


6161±1 


-2.0±0.3 


1000±200 


-10±2 






[O III] A5007 


7680±10 


1.2±0.2 


900±500 


4±1 






[O III] A5007 


7760±5 


0.8±0.2 


800±600 


8±2 


MP J2352-6154 


1.573±0.003 


C IV A1549 


3990±1 


0.7±0.1 


800±260 


80±16 






He II A1640 


4223±1 


0.40±0.04 


<500 


40±7 






C III] A1909 


4900±1 


0.30±0.04 


<500 


31±5 






[Ne IV] A2424 


6235±1 


0.40±0.05 


340±200 


30±10 






Mg II A2800 


7200±4 


0.3±0.1 


600±400 


30±10 



t Not considering the broad underlying component. 



